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I. INTRODUCTION 

^'. . . 

Oh' It is believed that strongly interacting matter under extreme conditions show a very rich phase structure. At 

^ high temperature and/or high density the hadrons overlap and loose their individuality; a new state of matter called 

(]J ' Quark Gluon Plasma (QGP) is formed ll|. It is well believed that such a state of matter existed in early universe, 

^S^ , a few seconds after the Big Bang, when the temperature was high enough for its formation. Deconfined quark 

matter should also exist in the core of neutron stars |2tJ| where the temperature is relatively low but density is high. 

So an understanding of the physics of strong matter in extreme temperatures and densities has cosmological and 

astrophysical significance. 

Such conditions of large temperatures and densities can be created in the laboratory by the collision of heavy ions 

'nJ" . at high energies. Relativistic heavy ion collision experiments include the Super Proton Synchroton (SPS) at CERN 

O^ ' in the past, at present the Relativistic Heavy Ion Collider (RHIC) at BNL and at Large Hadron Collider (LHC) in 

CERN and in the Facility for Antiproton and Ion Research (FAIR) at GSI. A plethora of data has been obtained 

t ' from RHIC, and a lot more is expected both from future runs there and from LHC, for strong matter at high energy 

^^ [ densities but nearly zero net baryon densities, which are relevant for cosmological considerations. FAIR, as well as 

future runs at RHIC, will be investigating matter at nonzero net baryon chemical potentials, and search for signatures 

of critical phenomena associated with a second order transition point. On the other hand, astrophysical laboratoris 

such as RXTE or EXOSAT may be expected to shed some light on the composition of matter inside a compact star. 

In this work we will concentrate on the strongly interacting matter at ultra-high densities and low temperature. As 

mentioned above, such a state of matter can exist near the core of a neutron star. Furthermore, inside the star the 
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d matter is charge neutral and at /3-equilibrium. Neutron stars are born at the aftermath of supernova explosions. The 

initial temperature of a new born neutron star can be as high as T ■^ lOOMel^. However, it cools down gradually 

mainly through neutrino emission and also through the emission of photons. For about one minute following its birth, 

the star stays in a special proto-neutron star state: hot, opaque to neutrinos, and larger than an ordinary neutron 

star (see, e.g., |5|,l6| and references therein). Later the star becomes transparent to neutrinos generated in its interior 

and transforms into an ordinary neutron star. The weak interaction responsible for the emission of these neutrinos 

eventually drive the stars to the state of /3- equilibrium and charge neutrality. The mass, radius and the composition 

of such a star depend on the nature of phase transition from hadron to quark matter [7|, |8[ . 

The possible central density of a compact star may be high enough for the usual neutron-proton matter to undergo 

a phase transition to some exotic forms of strongly interacting matter. Some of the suggested exotic forms of strongly 
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interacting matter are the hyperonic matter, the quark matter, the superconducting quark matter etc. Furthermore, 
there have been suggestions that the strange quark matter, containing ahuost equal numbers of u, d and s quarks, 
may be the ground state of strongly interacting matter Q. If such a conjecture is true, then there is a possibility 
of the existence of self-bound pure quark stars as well. In fact, the conversion of neutron star to strange star may 
really be a two step process [lO[. The first process involves the deconfinement of nuclear to two-fiavour quark matter; 
the second process deals with the conversion of excess down quarks to strange quarks resulting into a /3 equilibrated 
charge neutral strange quark matter. 

There are several mechanisms by which the conversion of strange quark may be triggered at the center of the star 
[111 Il2| . The dominant reaction mechanism for the strange quark production in quark matter is the non-leptonic 
weak interaction process |13j 

Ul + d -r^ U2 + S (1) 

Initially when the quark matter is formed Hd > fJ-s', and the above reaction converts excess d quarks to s quarks. But 
in order to produce chemical equilibrium the semileptonic interactions play an important role along with the above 
non-leptonic interactions 

d{s) ^ u + e^ + i7f, (2) 

u + er ^ d{s) + Ue (3) 

which implies the /3 equilibrium condition fJ.d=f^u+fJ'e+f^p] and ^s~lJ-d- 

Actually, the conserved charges in our system is baryon number and electric charge. Since we are assuming neutrinos 
to leave the system, lepton number is not conserved [ijj . Strange chemical potential fis is zero because strangeness is 
not conserved. So two of the four chemical potentials {nm fJ-d, y^si Me) are independent. In terms of baryon chemical 
potentiates) or quark chemical potential [iiq^iis/S) and charge chemical potential (/ig) these can be expressed as 

l^u= l^q+ gMg; fJ-d = fJ-q - -^^J■Q', fJ-s = ^J'q - ^fJ'Q', ^J■e = -fj-Q- 

Quantum Chromodynamics (QCD) is the theory of strong interaction and one should use this theory to look at 
the strongly interacting matter at high temperature and density. However, QCD is highly non-perturbative in the 
region of temperature and density that we are interested in. The most reliable way to handle the physics in this range 
is to use a lattice version of the SU{3) version of the colour gauge theory (lattice QCD). However, lattice QCD is 
numericaaly extremely costly and also not very reliable at finite baryon densities. So the most popular way is to use 
one of the effective models of QCD. Here we use one such model known as Polyakov - Numbu - Jona-Lasinio (PNJL) 
model [ll,[i3. 

It has been found that there is an instability in quark matter state with respect to diquark condensation, leading 
to a color superconducting state |17| . So at high density one expects a formation of diquark condensed phase. Here, 
in this paper, however, we have not considered this situation. 

In the present work we have studied the thermodynamics of /3 equilibrated 2-t-l flavour quark matter for non zero 
baryon chemical potenial. The study of /? equilibrated charge neutral quark matter has been carried out in NJL 
model earlier |18l. Il9j. Here the charge neutral trajectories are investigated in the T-^b plane. In the phase diagram 
of QCD we find the equilibrium phases consistent with charge neutrality and /3 equilibrium condition. 

Regarding the hydrodynamic evolution of matter, the s/ub trajectories are of importance. The behavior of entropy 
per baryon number in a plasma and in a hadron gas was analyzed within the framework of an extended Bag model by 
j2G |. The s/ub trajectories on QCD phase diagram was also studied using NJL [21[ and PNJL model [22]. In [23 the 
study has been done using the constraint of strange quark number density Ug = 0. Here, we study these trajectories 
considering (3 equilibrium conditions. The behavior of thermodynamic quantities like specific heat, compressibility 
etc. near the critical point is also studied for the system under consideration. 

This paper is organised as follows. In sec. II we discuss our model that we used. In the next section we calculate 
different thermodynamical properties and present our result and finally we conclude in sec. IV. 



II. FORMALISM 

The thermodynamic potential of 2+1 flavor PNJL model in finite quark chemical potential is [2j 
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Mf = m/ - 2gs<Tf + — cr/+icr/+2 
Electrons are considered as free non-interacting fcrmion gas and the corresponding thermodynamic potential is [1 
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where, fj,e is the electron chemical potential. An effective potential U{^,^,T) expressed in terms of the traced 
Polyakov loop $ = {TicL)/Nc and its (charge) conjugate $ = {Ti-cL^)/Nc is present in the PNJL Lagrangian. The 
Polyakov loop L is a matrix in color space explicitly given by 
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with j3 — \jT as the inverse temperature and A4 = iAP . The coupling between Polyakov loop and quarks is uniquely 
determined by the covariant derivative D^ in the PNJL Lagrangian. For simplicity, the temporal component of 
Euclidean gauge field A^ is treated as a constant in PNJL, and the Polyakov loop is reduced as 
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Corresponding to above expression, the trace of the Polyakov loop, $, and its conjugate, $, are treated as classical 
field variables in PNJL. 

The temperature dependent effective potential U{^,<^,T) is used to mimic pure-gauge Lattice QCD data, which 
should have exact Z{3) center symmetry. In this paper, we will use the potential U{^, $,T) proposed in [l^ , which 
takes the form 
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with 
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A precision fit of the coefficients Ui, hi is performed to reproduce the lattice data and Tq = 270MeV is adopted in our 
work. 

In Ref. [2^ a complete modification of effective potential (|9]) given by. 
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with J[$,$] = (27/247r2)(l - 6$$ + 4($3 + $3) „ 3($$)2) and k = 0.2. This J[$,$] is called Vandcrmonde 
determinant. The parameter set for the NJL part of our model is taken from Ref. j24| . 



III. RESULTS 



The thermodynamical potential fl, is extremised w.r.t the auxihary fields under the condition iJid=l^u-^l^e ', fJ-d—fJ-s- 
The equations of motions for the mean fields (7„, o^di <^si ^ a-nd 4> are determined through the coupled equations 
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This set of equations is then solved simultaneously for the fields cr,j, ad, (Jg ^ and $ as functions of quark chemical 
potential ^^ for different values of electron chemical potential ^e- 
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FIG. 1: (a) Quark constituent mass as a function of quark chemical potential at zero ^e ; (b)Quark constituent mass as a 
function of quark chemical potential at non zero /i^ 
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FIG. 2: (a) Quark number density as a function of quark chemical potential at zero ^e ; (b)Quark number density as a function 
of quark chemical potential at non zero /ie 



The number density of individual quarks and electrons are obtained by the relations 
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The total charge density nq — -^riu — hiT-d — k^s — ng. 

The pressure and entropy of the system is defined as P = — O and s = ^, respectively. The pressure of non-interacting 

quarks, gluons, and electrons is 
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FIG. 3: Total charge density as a function of quark chemical potential 
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FIG. 4: comparison of charge neutral trajectory in NJL and PNJL model at (a) jj,e — 10 ; (b) ^e=40 
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FIG. 5: (a) Pressure vs. inverse of number density ; (b) P/E vs. quark chemical potential at T =50 MeV. 



The specific heat C and compressibihty k can be evaluated using the definitions, 
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where, S is the entropy and n is the number density of the system respectively. In Fig(l) variation of quark constituent 
mass is shown as a function of ^q at T= 50 MeV. Both to„ and tti^, depending on the value of /ie, shows a jump around 
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FIG. 6: variation of compressibility k with quark chemical potential for (a) ^e=0, 10, 40 MeV and (b) 60 MeV; (b) Change 
in the minimum of P/E with fiq at T=50 MeV for different /le 



350 - 400 MeV region indicating a first order phase transition. On the other hand, such a jump will not be visible 
at higher temperature where transition becomes crossover. Moreover, since the equations for masses are coupled for 
non-strange and strange sector, the jump in m^ is smaller and is actually a manifestation of chiral transition in two 
flavor sector. This is indicated by the fact that ms remains again almost constant up to /ig=500 MeV and then drops 
sharply. 

In Fig 2(a) and (b) quark number densities arc plotted against fig at /ie=0 and 40 MeV respectively. For the first 
case, Uu = Tid and these become same as Us asymptotically for large fig. So for //e=0, total charge density uq of the 
system will be zero or the system will become charge neutral only for small fig or asymptotically at large fig as shown 
in Fig. 3. 

When fie is non zero, the fis and fid are greater than /j,„ due to /3 equilibrium. As we see in Fig. 2(b), at low fi region. 
Us is smaller than n„ as the higher mass of strange quark dominates the number density in this region. But at high 
fi this mass effect reduces and strange quark number density increases sharply and Ug > Uu due to higher chemical 
potential value. The electron number density is fixed for a fixed value of fig] and it is negligible compared to the Ug. 
At high fig region, there is a point where total charge is becomes zero due to the mutual cancellation of n„, n^, and 
ng, irrespective of n,.. Hence, Fig. 3 indicates the existence of charge neutral contour in the T-fig plane. For fi^ = 0, 
this contour is at ^^ = and at fig = oo. But for any non-zero fig, the contour exists for finite fig. These results 
can be reproduced for NJL model by putting (t>=(}>~\. In Fig. 4 the charge neutral trajectories for NJL model are 
compared with those of PNJL model, and it is seen that they are different in the hadronic phase, but in deconfined 
phase, the trajectories overlap. This is obvious because in PNJL model and is unity in deconfined phase. 

The system under investigation can be characterized by the behaviour of the EOS and shown in Figs. 5(a) and 
5(b). In Fig. 5(a) pressure increases sharply at lower values of n~^ or larger values of density. This variation of 
pressure with inverse number density indicates a repulsive behaviour of the interaction in the many body system at 
large densities (large fLg) or short distances and an attractive nature at larger distance or lower densities [26|, |27[. In 
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FIG. 7: Variation of specific heat scaled by its Stefan Boltzmann value, with fiq at 
MeV for different electron chemical potential 



a) T=50 MeV; (b) T=90 MeV; (c) T=140 



Fig. 5(b), P/E is plotted against fj,q at T=50 MeV for various fie- The ratio P/E starts to decrease around fj,q= 250 
MeV and reaches a minimum around /iq=350 MeV. This change is indicative of the softening of the EOS in first order 
phase transition. This in turn corresponds to the reduction in isoentropic speed of sound. In a 1st order transition, 
in the presence of one conserved charge, pressure is constant in the mixed phase and speed of sound (dP/dE) is zero. 
On the other hand, for two conserved charges, pressure is not constant any more, rather its variation becomes slower 
[28| resulting in a smaller but non-zero speed of sound. The second smaller peak at around iig~ 450 MeV in Fig. 6(a) 
is due to the mass change in the strange sector. 

The signal of a first order phase transition and the small shifts in the transition point for various fie can be seen in 
the variation of compressibility, P/E, and specific heat plotted in Figs. 6(a), 6(b), 6(c), 7 (a, b, c) respectively. Both, 
compressibility as well as specific heat are second derivative of fl and represent the fluctuations near the transition 
point ,26(]. Discontinuity in compressibility as well as specific heat at T=50 MeV indicates first order phase transition 
for fie=0, 10, and 40 MeV. Various values of fig for different fie , in this region, then correspond to the values of quark 
chemical potential in the mixed phase. For /ie=60 MeV, transition is a crossover at this temperature. 

In Fig. 8; the contour of constant baryon number density, scaled by the normal nuclear matter density (no ~ 
0.15fm~^) is shown along with the phase diagram. The charge neutral trajectories for each non zero ^e are also 
shown. At ij.q ^ 450 MeV, the density is about 5-10 times the normal nuclear matter density. Such densities are 
expected to occur in the neutron star core. Only for /ig ^ 40 MeV, (Fig 8c, 8d) the charge neutral trajectories lie 
in the regions of such large density and low temperature. This result gives an estimation of the value of electron 
chemical potential consistent with the properties of neutron star. 

The net strangeness fraction {ns/ns) along with n^/no is shown in Fig. 9. This figure along Fig. 8 indicates that 
at the core of the star the strangeness fraction will much larger and one can expect the formation of quark matter 
with almost equal number of u, d and s quarks. Similar results have also been found in other model studies [29|. 
Moreover, the intersection of lines of constant baryon density and strangeness fraction indicates the possibility of 
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evolution of the system to higher strangeness fraction at a constant density. In hydrodynamic simulation the entropy 
per baryon number density ratio {s/ub) is an interesting quantity. This is because the adiabatic hydrodynamic 
expansion conserves s/ub along the time evolution direction. This prescription is appropriate for discussing the 
thermodynamics of matter created in heavy ion collision, because after equilibration the fireball will expand along 
the line of constant s/ub- On the other hand, the scenario inside neutron stars will be much more complicated, as 
discussed later, due to much larger time scale involved. 

By studying s/ub one can check that if there is a converging tendency of the adiabats towards the critical points, 
as claimed in J3C| . We analyze here how the quark chemical potential changes along the constant s/ub line, under 
flavor equilibrium constraint. In Fig. 10 (a, b, c) the isentropic trajectories are shown for 2+1 flavour quark matter 
for various /ig. As we see the nature of the trajectories is quite similar to those of [21| and [2j|- At T ^ 0, s — )• by 
the third law of thermodynamics; in order to keep s/ub constant, ub should go to zero. This condition is satisfied 
when /i — i- Myac of the theory [2ll |. and we see the curves have a tendency to converge at that point. On the other 
hand, in Fig. 10(d), the isentropic trajectories are shown for 71^=0 and /ie=0. The constraint rig =0 is necessary to 
study the matter to be created in CBM experiment, it is basically a 2 -flavor matter. However, there is no significant 
difference between the trajectories of Fig. 10(a) and 10(d). We can expect that the experimental results should agree 
with our model prediction for 2+1 flavour case. 

In the light of above results, let us now discuss phase transition inside neutron stars. Since the temperature inside 
a new born neutron star drops very quickly, one may assume it to be a system of low temperature nucleonic matter 
system which may also be populated by hyperons and strange baryons due to high density near the core. Let us 
also assume it to be /3 equilibrated and charge neutral. Now due to some reason {e.g. sudden spin down) conversion 
may start and this nucleonic matter will get converted (or deconfined) to predominantly two flavour quark matter 
within strong interaction scale. This conversion may start at the centre and the conversion front moving outward will 
convert most or all of the star. Each point inside the star is expected to lie on isentropic trajectory. This system of 
predominantly 2 flavour quark matter then will get converted to strange quark matter through weak interactions and 
finally a j3 equilibrated charge neutral strange quark matter will be produced. This 2+1 flavor strange quark matter 
is more stable as we see from Fig 11. Here the pressure normalised by its Stefan-Boltzmann value, is compared for 
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2-flavor and 3-flavor quark matter. The bump at small ^.q arises due to the eleetronic eontribution to the pressure. 
Since P=-i7, the free energy of the 2+1 flavor system is smaller than that of 2 flavor system, hence more stable. The 
strangeness production occurs mainly through non-leptonic decay [IJ], the system is expected to lie on a constant 
density line and move towards the point with highest strangeness possible at that density. Finally the semi-leptonic 
processes will take over and system will then evolve along a j3 equilibrated charge neutral contour. 



IV. CONCLUSION 



In this paper we have studied the strongly interacting matter at moderate temperature and high density under 
beta equilibrium. For this purpose we have used 2+1 flavor PNJL model in the mean field approximation. The 
thermodynamic potential has been calculated and it has been minimised with respect to the fields to obtain the field 
values at a particular temperature and density. Once the field values have been obtained these values are then used 
to calculate the constituent quark masses. We have looked at the masses for both finite and zero electron chemical 
potential and at T = bOMeV. At this value of temperature phase transition from hadron to quark matter is first order 
for both zero and finite electron chemical potential. The quark number densities have been studied for same values 
T and /ig and those are found to increase steeply with density. It is well known that inside a neutron star the matter 
is charge neutral. So we have studied the charge neutral contour of beta equilibrated matter in the T — fis plane and 
compared our results with the NJL model. We have found that the nature of the contour depends strongly on the 
value of the /ig. The equation of state has been studied through the calculation compressibility and P/E. The minima 
of the P/E is found to be quite sensitive to the electron chemical potential. The compressibility also depend on the 
electron chemical potential strongly. So one can conclude at a moderate temperature the net electron density governs 
the eos and also the order of the phase transition which can be confirmed by the study of specific heat. We have 
obtained the charge neutral contour separately and found the equilibrium phases consistent with charge neutrality 
and /3 equilibrium conditions. The value of electron chemical potential is taken care of accordingly. We have also 
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FIG. 10: The isentropic trajectories along with phase diagram at (a) He—O ; (b) /^e=40; (c) /ie=60; (d) 7X3=0, fJ,e=0', The curves 
correspond to (from left) s/7is=300,100,30,10,5,3.5 
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FIG. 11: Comparison of pressure for 2-flavour and 3-flavour quark matter for fie—O 



obtained the isentropic trajectories along which the system is expected to evolve. 

In this work we have studied the phase diagram keeping the electron chemical potential as a fixed quantity. The 
natural extension of the work is to obtain the phase diagram by varying both the quark and electron chemical 
potential. Moreover, since /3 equilibrium and charge neutrality are two major conditions for constructing a neutron 
stars, present results may be applied to study the neutron stars configurations. Finally, the work has been done without 
incorporating di-quark condensate. It would be interesting to include this in our model and study its consequences 
for neutron stars. 
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